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Abstract. I describe the logical basis of the method that I developed in 1962 
and 1963 to define a quantum operator corresponding to the observable particle 
| number of a quantized free scalar field in a spatially-flat isotropically expanding 

■ (and/or contracting) universe. This work also showed for the first time that particles 

. were created from the vacuum by the curved space-time of an expanding spatially-flat 

O { FLRW universe. The same process is responsible for creating the nearly scale-invariant 



spectrum of quantized perturbations of the inflaton scalar field during the inflationary 



$H , stage of the expansion of the universe. I explain how the method that I used to 

obtain the observable particle number operator involved adiabatic invariance of the 
particle number (hence, the name adiabatic regularization) and the quantum theory of 
measurement of particle number in an expanding universe. I also show how I was led 
in a surprising way, to the discovery in 1964 that there would be no particle creation 
\^ , by these spatially-flat FLRW universes for free fields of any integer or half-integer spin 

satisfying field equations that are invariant under conformal transformations of the 
metric. The methods I used to define adiabatic regularization for particle number, 
were based on generally-covariant concepts like adiabatic invariance and measurement 
that were fundamental and determined results that were unique to each given adiabatic 
order. 



1. Introduction 

It is my great pleasure to write this article for a special issue of J. Phys. A: Math. 
Gen. in honor of Professor Stuart Dowker. I am particularly grateful to him for his 
outstanding contributions to quantum field theory in curved space-time, from which I 
first learned about the beautiful method of zeta-function regularization in curved space- 
time PEISIIH]. 

It seems particularly appropriate for me to describe here the regularization 
technique that I developed in 1962-1964, while working on my Ph.D. thesis [6], to define 
an operator corresponding to the observable particle number of a free scalar field in a 
spatially-flat isotropically expanding (and/or contracting) universe. This construction 
appears in detail in |6], but until now I have not discussed it fully in the literature. I 
hope this discussion will clarify the basis and uniqueness of the construction, as well 
as its intimate relation to particle creation by the gravitational field of the expanding 
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universe. As you will see, the method is not limited to weak gravitational fields or to 
any particular form of the change in the scale factor of the universe, a(t), beyond the 
continuity of a(t) and of a sufficient number of its derivatives. 

The need for regularization of the particle number observable became clear to me 
only after I had considered a quantized field in a smoothly expanding spatially-flat 
FLRW universe that started from a Minkowski space at early-times and ended in a 
Minkowski space at late-times. There was no question how to define the vacuum state 
and particle number of a free quantized field in Minkowski space, and there was no doubt 
how to propagate the field into the curved space-time of the expanding universe using 
the generally-covariant equation for the free field obtained by the methods of general 
relativity. Therefore, I was able to show convincingly that particles were definitely 
created from the early-time Minkowski vacuum state by showing that after the expansion 
of the universe there were particles present in the late-time Minkowski space. The 
gravitational field of the expanding universe was regarded as a classical field while 
the quantized field was treated using quantum field theory propagated into the curved 
space-time of the expanding universe by the generally-covariant free field equation. 
Perturbations, such as that of the gravitational field (and other classical fields) satisfy 
linearized field equations, so the same methods also apply to such perturbations, which 
can be treated as quantized fields that are also created by the expansion of the universe. 
For a statically bounded (or asymptotically static) smooth expansion of the spatially- 
flat FLRW universe, I found that no regularization was needed, apart from the usual 
Minkowski-space vacuum subtractions at the beginning and end of the expansion, to 
prove that in a co-moving volume (1) the particle number was an adiabatic invariant, and 
(2) that a smooth generic statically bounded expansion of the universe, starting from the 
initial Minkowski space vacuum state, would create particles as a result of the expansion 
of the universe, giving a non-zero expectation value of the particle number operator in 
the late-time Minkowski space. The total number density of created particles, summed 
over all modes of the quantized field, was finite. 

Let me digress for a moment to point out that this particle creation is the same 
process that is responsible for the creation of the nearly scale-invariant power spectrum 
of primordial inflaton perturbations created by the inflationary stage of the inflationary 
model of the early universe. A nice proof was given in [7|. We proved this by 
joining an asymptotically static smooth expansion of a spatially-flat FLRW universe 
that approached a Minkowski space at early times to an arbitrarily long period of 
exponential inflation that was finally joined to another asymptotically static smooth 
expansion of the spatially-flat FLRW universe that approached another Minkowski 
space at late times. The joining of the metrics of these three stages of the expansion 
were continuous to second-order in time derivatives of the expansion scale (as that 
was sufficient to have no UV divergences in the created number of particles in any co- 
moving spatial volume). Furthermore, the solutions of the wave equation were required 
to have appropriate continuity conditions on their time derivatives. (This problem is 
mathematically analogous to solving a time-independent quantum mechanical scattering 
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problem in one spatial dimension.) The result was calculated analytically and then 
evaluated to many digits of accuracy. Our result confirmed that if we started from 
the Minkowski vacuum state of the inflaton perturbation field at early times, then 
the state evolved after a relatively small number of e-foldings into the Bunch-Davies 
vacuum state of the inflationary universe (except for wavelengths that were so large 
that they were already outside the inflationary Hubble radius when inflation started and 
remained outside when inflation ended). The result was that the spectrum of inflaton 
perturbations by the end of 50 or more e-foldings of inflation was scale invariant for the 
range of wavelengths one could observe in the present universe. Because the expansion 
of the universe was asymptotically static in the past and the future, there were no UV 
infinities that survived to late times. There were also no infrared (IR) divergences. 
The main point I am making here is that the creation of particles (or quanta of a 
field) by the expansion of the universe is responsible for producing the spectrum of 
inflaton perturbations in the standard inflationary scenario. Thus, the agreement of 
current observations of the CMB power spectrum with the inflationary prediction for the 
primordial spectrum of inflaton perturbations can be interpreted as evidence in favor of 
the particle creation process we are discussing in the present paper. 

Furthermore, if inflation did begin from some kind of initial state that was not the 
vacuum, then we may be able to observe some features of the initial state prior to the 
beginning of inflation, even if inflation lasted for as many as 60 e-foldings or more. This 
goes back to another topic that appeared in my Ph.D. thesis and early papers, namely, 
the stimulated creation of bosons by the expansion of the universe [6j [161 E] Recently, it 
was shown that this effect may show up in observations of non-Gaussianities that may 
appear in the bispectrum of CMB anisotropics [Hi [19] , as well as in the trispectrum 

Let me now return to our discussion of the method of regularization that I developed 
in 1963 and 1964 to obtain a suitable particle number operator for a free quantized field 
in an expanding universe. Encouraged by the results I had obtained for a statically 
bounded expansion of the universe, I went on to consider universes that were not 
statically bounded. After all, observation had shown that we live in a stage during 
which the universe is expanding. However, if I tried to use this same time- dependent 
particle number operator to define the number of particles present at a time t during 
which the universe was undergoing a generic expansion with non-zero time-derivatives 
of the scale factor a(t), then I found that the average number of particles (in a given 
co-moving volume) when summed over all modes would be infinite, as the result of a 
UV divergence proportional to the square of a(t)/a(t) and to a{t)/a{t). 

The particle number observable of a quantized field in the curved space-time of 
an expanding universe has certain physically-motivated requirements and limitations 
that guided my reasoning. The first physically-motivated assumption, underlying the 
method of regularization, was that the particle number of a free (non-interacting) 
quantized field in the limit of an infinitely slow and smooth expansion of the universe 
should not change with time. This assumption that the particle number is an adiabatic 
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invariant is one fundamental requirement of this regularization method, which therefore 
is known as the method of adiabatic regularization. (It is clear that requiring adiabatic 
invariance does not imply that the particle number operator we obtain will be valid 
only for scale factors a(t) that change slowly relative to parameters associated with 
the free quantized field and the gravitational model under consideration.) A second 
physically-motivated assumption fundamental to the regularization method was that 
the regularized particle number observable must have the non-negative integers as its 
eigenvalues. (A measurement of the particle number present at a given time may be 
carried out over an interval of time, just as a measurement of the momentum that a 
particle had at a given time may require an interval of time to complete.) A third 
requirement was that during any period when a(t) was constant, the particle number 
observable would reduce to the particle number observable that I had used in the 
statically-bounded expansion, and which I already had proved was an adiabatic invariant 
for such statically-bounded expansions. 

In later Sections of this paper, I will describe how I found a particle number operator 
that satisfied these requirements and gave a finite total number of created particles per 
unit co-moving volume when summed over all modes of the field at any time during the 
expansion of the universe, even if the expansion was not statically bounded. In defining 
the particle number operator, I made the minimal number of adiabatic subtractions that 
were necessary to obtain a finite total number of created particles per unit co-moving 
volume, while satisfying the requirements just described. (Some years later, in 1974, 
S.A. Fulling and I used this method of adiabatic regularization [HI [9] to renormalize the 
energy-momentum tensor, which required one more term in the adiabatic series for its 
regularization, as compared with the adiabatic regularization of the particle number.) 

In the course of obtaining these adiabatic subtraction terms, I found that there was 
a special equation of motion of the field for which there was no creation of real particles 
for any given smooth expansion of the universe. (By saying real particles here, I exclude 
virtual pairs of particles that would contribute to vacuum energy and pressure during the 
expansion of the universe, but would not be present at late times in a statically bounded 
expansion of the universe.) I showed that particles satisfying this special equation of 
motion had a unique well-defined vacuum state throughout any smooth expansion of 
the universe. By comparing that special equation of motion with the actual equation 
of motion of the particles under consideration, I was able to obtain expressions for the 
subtraction terms that satisfy the requirements stated in the previous paragraph. 

I had completed this work in 1963 for the scalar field and during 1964 was in the 
process of extending these results to spin-1/2 fields, when I happened to see a 1964 article 
of Roger Penrose [10] in which he wrote down the particular wave equation for a scalar 
field that was invariant under conformal transformations of the metric. I was startled 
to see that this equation was exactly the same as the special equation for the quantized 
scalar field that I had earlier found gave rise to no creation of particles and had a unique 
well-defined vacuum state throughout any smooth expansion of the spatially-flat FLRW 
universe. 
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The remarkable identity of these two equations, led me to find a straightforward 
proof that there was no particle creation for free fields of any integer or half-integer 
spin, that satisfied the conformally invariant field equations in the spatially-flat FLRW 
universe. This proof was completed in 1964 and was included in my Ph.D. thesis. 
Among the conformally-invariant field equations are the massless Dirac equation and 
the Maxwell spin-1 photon field equation. However, one should note that the spin-2 
linearized graviton equation obtained from the Einstein gravitational field equations is 
not the same as the conformally-invariant spin-2 field equation. 

In the following Sections, I will try to elucidate these introductory remarks relating 
to the discovery of the creation of particles from the vacuum in an expanding universe 
and the related development of the method of adiabatic regularization of the particle 
number observable in such a universe. 

2. A Brief Orientation 

In 1962, when I began my thesis work at Harvard under the direction of Sidney Coleman, 
my objective was to study the possibility that particles of quantized fields in the curved 
space-time of an isotropically expanding (FLRW) universe would be created from the 
vacuum as a result of the expansion (or contraction) of the universe. At the time, I was 
not aware of some previous work by Schrodinger related to this topic [TT| [12]. A brief 
discussion of his work is given near the end of Section [91 

In curved space-time, I soon realized that methods like normal ordering that set 
the energy of the vacuum state to were not acceptable because the energy density 
of the quantized field must appear in the Einstein gravitational field equations, and if 
particle creation by the expanding universe does take place, then the energy density of 
the created particles would not be expected to remain as time evolved, even if one 
started with an initial vacuum state of energy. To understand what must occur if 
there is particle creation, let us use the Heisenberg picture, in which the state vector 
does not evolve with time, but the quantized field operators do evolve. 

In the Heisenberg picture, suppose that the state vector describing the free 
quantized field in the expanding universe contains no particles of the field at early times. 
How may the quantized field evolve with time in such a way that the same state vector 
at late times will be found to contain a non-zero number of particles of the field? Such 
an evolution would require that the annihilation operators of the early-time quantized 
field evolve into linear combinations of annihilation and creation operators for quanta of 
the late-time field. (Otherwise, the annihilation operators of the late-time field would 
annihilate the unchanging state vector, in which case there would be no particles present 
at late times and the initial vacuum state would also be the final vacuum state). So if 
particle creation does take place, one must have time evolution under the Heisenberg 
equation of motion, such that as the quantized field, say <p evolves in time, the late time 
annihilation operator is a linear combination of the initial annihilation and creation 
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operators and A_t\ 

B % =a k A % + p k Aj. 

This transformation corresponds to the creation of particle-antiparticle pairs of total 
momentum from the vacuum, as required by conservation of momentum. This linear 
transformation is an example of a Bogolubov transformation, as was first explicitly 
pointed out in the context of gravitationally-induced particle creation by [16J. Analogous 
Bogolubov transformations that involve a sum over different momenta, as in 

B k = J^Kfc' A %> + Pk,k> A_^]. 
P 

occur in other situations that involve gravitationally-induced particle creation. 
References may be found, for example, in [13] and [T4] . 



3. Free Scalar Field in Expanding Universe 

Let us now consider the case of a free quantized scalar field in a smoothly changing 
spatially flat FLRW universe with line element 

ds 2 = dt 2 - a 2 (t)(dx 2 + dy 2 + dz 2 ), (1) 

and equation of motion 

(□ + m 2 )0 = O, (2) 

where □ = /"V^Vy, with V/ denoting the covariant derivative. For now we are 
considering the so-called minimally- coupled scalar field having mass m. More generally, 
one could include a coupling to the Ricci scalar curvature, R = g^ u R^ U) namely, 
(□ + m 2 + £R)(j) = 0, where £ is a dimensionless constant. We work in units with 
H and c each equal to 1. 

Writing ([2]) with the metric of (pQ), one finds that 

a~ 3 d t (a 3 d t <j)) - a" 2 ^ 9f0 + m 2 = 0. (3) 

i 

It is convenient to impose periodic boundary conditions in a cube having sides of 
coordinate length L and coordinate volume V = L 3 . As in Minkowski spacetime, this is 
a mathematical device, with L taken to infinity after physical quantities are calculated. 
Then we can expand the field operator <fi in the form 

k 

where 

f s = {2Va{tf)- l ^ s h k {t). (5) 
Here k % = 2nn l /L with n l an integer, k = \k\, and h k (t) satisfies 

^h k + ^h k {t) + m 2 h k (t) - - A (-) h k (t) - \-h k {t) = 0. (6) 



dt 2 a 2 K ' w 4 V a / " w 2 a 
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We infer that 




(/*#) = o. 



(7) 



where the quantities on the left-hand-sides of <^ denote scalar products that are 
conserved by virtue of ([3]). These scalar products involve integration over a spatial 
volume at time t in the FLRW universe. (On the conserved scalar product and its 
consequences, see [H] [Section 2.2, (2.48)].) The right-hand-sides of (JTj) are necessary to 
insure that if the scale factor a(t) were to smoothly approach a constant, for example, at 
late times, then the scalar product would agree with that of a free field in the late-time 
Minkowski space. (Consistency requires © to hold even if a(t) were not to actually 
approach a constant in the distant future or for any particular time-interval.) The 
constant creation and annihilation operators in ((1]) satisfy the commutation relations, 



For the relation of these to spin and statistics in an expanding universe, see [14] . 

In the next Section, I introduce time-dependent creation and annihilation operators 
as I did in [6j[for bosons and fermions], [15] [for bosons and fermions], [16j[for bosons], 
and [17] [for fermions]. I also discuss the adiabatic invariance of particle number, and in 
a later Section, the role of measurement of particle number in determining the adiabatic 
subtraction terms. The creation of particles by the expanding universe in the context 
of quantum field theory was discovered in conjunction with the previous ideas. The 
connection of particle creation with conformal invariance was also discovered in the 
course of developing these ideas in my Ph.D thesis. These topics will arise naturally in 
the following Sections. 

4. The origin of the method of adiabatic regularization 

In order to help the reader understand the original physical basis for adiabatic 
regularization, it is necessary for me to explain some crucial features that appear in 
detail in my Ph.D. thesis [6]. In the early papers that I published on the results in [6], 
my emphasis was on particle creation by the gravitational field. The method I developed 
to regularize the time-dependent particle number during the expansion of the universe 
was only briefly outlined. 

Here, I will emphasize the basic ideas and go into sufficient detail to show how the 
regularization was based on generally-covariant concepts, such as adiabatic invariance 
defined in terms of proper (cosmic) time t along the preferred geodesies in an FLRW 
universe, and measurement of the number of particles present in a co-moving volume 
at a given cosmic time. If such a measurement is made over a period of time when 
particles are being created by the expansion of the universe, then there is an irreducible 
uncertainty in the result of the measurement. This uncertainty allows sufficient latitude 
in the definition of the particle number operator to allow us to regularize the original 
number operator so that it remains hermitian, has non-negative integer eigenvalues, and 



[Ak,^]=0, =0, and A h A\, 




(8) 
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has no UV infinities when summed over all modes of the field. The uniqueness of the 
resulting particle number observable, defined up to a given adiabatic order, is briefly 
explained, as is the minimum adiabatic order (in this case second adiabatic order) 
required to yield a particle number operator that has no UV divergence during the time 
when particles are being created. If one considers a statically bounded expansion of the 
universe, then the total number of particles created as a result of the expansion is the 
same as I had obtained earlier without regularization. 

These features are well illustrated in the spatially-flat FLRW universe. Unless 
otherwise noted, the term FLRW universe in this paper refers to the class of spatially- 
flat FLRW universes with a(t) smoothly changing. 

There are several basic assumptions that guided my reasoning. The first is that 
the particle number in a smoothly expanding universe should be an adiabatic invariant. 
This means that in the limit of an infinitely slow and smooth expansion (the so-called 
adiabatic limit), the particle number should remain constant, in the sense that the 
change in the particle number should approach zero in the adiabatic limit even if the 
change in the scale factor that occurs over a very long time, i.e., the ratio a(t2)/a(ti), 
is very large. The adiabatic approximation of Liouville for the harmonic oscillator with 
a slowly changing frequency gives the following approximation to the solution of (|6]): 

h k (t) ~ (u; fc (t))- 1/2 exp(±z / u k {t')dt'), (9) 



where 

u k (t) = v/(fc/a(t)) 2 + m 2 . (10) 

This approximation is good in the limit that all time-derivatives of a(t) smoothly 
approach o|j| The two solutions in ([9]) are linearly independent. Then the general 
solution of the second-order ordinary differential equation OH]) can be written in the 
adiabatic approximation clS db linear combination 



hk(t) ~ a k (u k (t)) 1/2 exp(-i J u k (t')dt') 

+ P k {u k {t))~ 1 / 2 exp{+t j u k (t')dt'), (11) 

where a k and (3 k are complex constants that must satisfy 

K| 2 -|/3 fc | 2 = l (12) 

because of ((Tj). In general, a k and (3 k will have a time-dependence that is determined 
by the solution of the ordinary second order differential equation ([6]) with h k (t) written 
as the sum of two linearly independent solutions in the form 

h k (t) = a k (t) {u k (t))- 1/2 exp(-i / w*(0 dt') 



| The Liouville approximation is closely related to the JWKB approximation developed later to 
approximate solutions of the Schrodinger equation for slowly changing scattering potentials, including 
the possibility of tunneling. 
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+ /3 k (t)(uj k (t))- 1 / 2 exp(+t J u k (t')dt'), (13) 

with 

\a k (t)\*-\p k (t)\ 2 = l. (14) 



We will assume that a(t) and all its time-derivatives are smooth and well-defined. 

Now the requirement that the particle number in each mode of the field <f>{x, t) 
should be an adiabatic invariant is easily met by factoring out the positive and negative 
frequency adiabatic approximations that appear in ( TTTT) and including the remaining 
time-dependence of the exact solutions of ([6]) in the definition of time-dependent creation 
and annihilation operators. Thus, we rewrite (jlj) in the form 

k 

where 

g R (x,t) = V- 1 / 2 e^(2a(t)- i u k (t))- 1 / 2 exp(-t J* u k (t') dt>). (16) 

and 

a^t) = a k (t) A^ k + (3 k (t)* A^. (17) 

As a consequence of (Tl2l and other conserved quantities (as discussed in [Ej), it follows 
that 

= %m (18) 

Thus, the al(t) and aAt) satisfy the commutation relations of creation and annihilation 
operators at all times t. 

Can these operators a-t(t) and a k {t) be interpreted as annihilation and creation 
operators for real particles present at time t? This interpretation is well-supported 
and consistent at early and late times when one considers a statically bounded (or 
asymptotically static) expansion of the universe for which a(t) smoothly approaches a 
constant value a\ at early times and a constant value a-i at late times. When I began my 
Ph.D. thesis work in 1962 this was the appropriate case to look at in order to decide if 
physical particles were created by the expansion of the universe. The reason for having 
the geometry approach Minkowski spaces at early and late times was that one knew 
unambiguously how to define real particles in Minkowski-space quantum field theory. 
Furthermore, that interpretation agreed with experimental results in Minkowski space. 
I was able to prove by 1963 that real particles were indeed created as a result of any 
typical statically bounded smooth expansion (or contraction) of the universe. I also 
showed that the number of particles was an adiabatic invariant, in the sense that the 
number of particles present at late times was the same as the number that had been 
present at early times in the limit of an infinitely slow and smooth expansion of the 
universe P SSI UM [12] . 



[a%(t),ap(t)] = 0, 



0, and 



a fc(*)>4^ 
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In a later Section, we will take up the interpretation of these operators a>r(t) and 
during a time when a(t) is smoothly changing. As we will see, the particle number 
in generic cases involves a UV infinity when summed over all modes of the field. This 
UV infinity goes away at early and late times for the class of smooth statically bounded 
expansions of the universe, so no regularization was required to show that particles were 
created by a statically bounded expansion. However, after showing with no ambiguity 
that particles are created in such statically bounded expansions, it became necessary to 
consider the meaning of particles during a time like today, when the universe is actually 
expanding. 

In my thesis, I developed a method based on an asymptotic series expansion about 
the adiabatic solution of the field equation to deal with the UV infinities that appear 
in expectation values of the operator at(t)a^(t). The physical basis for the method 
was a careful consideration of the measurement of the particle number at a given time 
when a(t) is smoothly changing as a function of time. This was the original basis for 
adiabatic regularization. In the context of particle number its development was guided 
by the natural requirement that the particle number for the free field should be an 
adiabatic invariant and that the measured particle number should be an integer. 

We have seen in this Section how the adiabatic condition arose in a natural way 
when considering particles and their creation by the expanding universe; and we will 
see in a later Section the role played by the adiabatic condition in the construction of 
a suitable number operator corresponding to the particle number measured during the 
expansion (or any smooth change) of the universe. I am often asked why adiabatic 
regularization gives covariant results. The reason is that the concepts of adiabatic 
invariance, and the requirements we discuss later that must be satisfied in a measurement 
of particle number in an FLRW universe, are independent of coordinate systems, being 
based on proper time t along the preferred congruence of geodesies assumed to exist 
in the FLRW universes. To make this clear to the reader, I will include below a fairly 
detailed discussion of the considerations concerning measurement of the particle number 
in an expanding universe. These considerations appeared in detail in my thesis, but were 
only briefly sketched in my early publications on my thesis ([IS]). Therefore, I think 
they will be novel and of value to the reader. 

Before discussing adiabatic regularization of UV infinities in the expectation value 
of the particle number summed over all modes, we introduce some background and 
discuss an important relationship discovered in my thesis that involves particle creation 
and a preferred vacuum state having a symmetry additional to the usual symmetries of 
the FLRW universes. 
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5. Scalar field equation satisfied exactly by a general adiabatic-type 
solution 

The differential equation having exact solutions of the form 
h k (t) = a k (W k (t))- 1 / 2 exp(-t j W k (t')dt') 

+ P k (W k (t))- 1 / 2 exp(+t J'w^dt'), (19) 

with both a k and (3 k independent of time, is 

h(t)o + {W k {tf - W^Kd'/dt^W^ty 1 / 2 ]} h k (t) = 0. (20) 

We will assume that the function W k (t) is continuous and has at least two continuous 
t ime- deri vat i ves . 

For the minimally-coupled scalar field of ([2]), we can rewrite (El) in the following 
form [6J[p.28]: 

h k (t) + {co k (t) 2 - u k {tfl 2 [{d 2 /dt 2 )u k {t)- l / 2 }}h k {t) = 2u k {t)S{k,t)h k {t), (21) 

where u k (t) given in (flOl . Then, the function S(k,t) for this minimally-coupled scalar 
field is defined by 

2u k (t)S(k,t) = dik^^d^/ait)) 2 + C 2 (k,t)(a:(t)/a(t)), (22) 

with 

fc 4 + 3m 2 a(t) 2 k 2 + (3/4)m 4 q(t) 4 
Cl{k > t] ~ (k 2 + m^aitf) 2 ( j 

and 

k 2 + (3/2)m 2 a(t) 2 

° 2 ^- V + mMt) 2 { ) 

As an aside, one can show that 2u k (t)S(k,t) vanishes in the non-relativistic limit 
that [m 2 a 2 /k 2 ) 3> 1 for the particular solution, a(t) oc t 2 ^ 3 , of the classical Einstein 
equations for a spatially-flat FLRW universe dominated by non-relativistic particles. 
In addition, one can show that 2u k (t)S(k,t) vanishes in the relativistic limit that 
(m 2 a 2 /k 2 ) 1 for the particular solution, a{t) oc t 1 ^ 2 , of the classical Einstein equations 
for a spatially-flat FLRW universe dominated by relativistic particles. This was first 
pointed out in [6j [15], [16] and means that in those universes there is no creation of 
these particles of the dominant type. This is because the adiabatic solution is the exact 
solution of the minimally-coupled free field equation in these two cases, which means 
that the \(3 k \ = for them, so there is no particle creation of the dominant type of 
particle. 

Later we will make use of f l2~T|) in arriving at the adiabatic subtractions that 
regularize the UV divergences that appear in the total particle number summed over all 
modes when a(t) is not constant. 

But first we use ( 12T|) to show that there is a particular non-minimally coupled 
scalar field equation having a definite preferred vacuum state over the entire period 
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during which a(t) is smoothly changing. I found this result in 1963 and it led me to the 
discovery of a surprising connection between particle creation and a symmetry that I 
only became aware of a year later in 1964. 



6. Scalar field equation giving no creation of massless particles by any 
statically bounded smooth expansion of the spatially- fiat FLRW universe 



In the case when m = in (Izlj) . we find from (|22j) - (j24|) that 

h(t) + {oo k {tf - u k (ty/ 2 [(dy deputy 1 / 2 ]} h k (t) = (i/6)R(t)h k (t), (25) 

where R(t) is the Ricci scalar curvature in the spatially-flat FLRW universe, 

R(t) = 6 [(d(t) I 'a{t)f + (a(t)/a{t))] . (26) 

From this we immediately deduce that if we were to start from the massless scalar 
field equation similar to the minimally coupled one in (J2J), but with an additional term 
involving (l/6)i?(t), namely 

(□ + (1/6)^)0 = 0, (27) 

which in the spatially-flat FLRW metric is 

a~ 3 d t (a 3 d t (f)) - or 2 d2 i<t> + (1/6 W = °> ( 28 ) 

i 

then the additional (1/6) R on the left-hand side of (I28I) would cancel the same term 
that appeared on the right-hand-side of ( 1251) when we had started from the minimally- 
coupled massless scalar field. Therefore, if we were to start from (1271) instead of (j2J), 
then the same steps that had led to ( 1261) would lead us to the result that 

h(t) + {Mt) 2 - ^(0 1/2 [(d7^>*(r 1/2 ]} h(t) = 0. (29) 

But this is the same as (I20p . the equation that earlier we noted has as its exact general 
solution (119j) . which in the present case takes the form 

h k (t) = Oi k (w fc (i))~ 1/2 exp(-z J u k (t') dt') 

+ /3 fc (u; fc (t))- 1/2 exp(+2 J u k {t')dt'), (30) 

with a k and j3 k any chosen complex numbers independent of t and satisfying ( 1T21) . 

It is now easy to go over the discussion in Section [31 but with (T5]) replaced by (127]) 
and ((6]) by ( 1291) . to see that if we choose a k = 1 and (3 k = in ( 1301) . then a k (t) = A k 
and a^ k (t) = A\ for all t, which implies that there is no creation of real particles in any 
smooth statically-bounded expansion. Consistency then implies that there is no creation 
of real particles satisfying this free, massless, scalar field equation, (127]) . for any smooth 
non-singular expansion a(t), even if there are no periods when a(t) is not changing with 
time. Therefore, for the field satisfying ( 1271) there exists a preferred vacuum state having 
no real particles at any time during the expansion. This preferred vacuum state is the 
state annihilated by the operators A k for all k. 
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The presence of interactions with other particle fields or of self-interactions would, 
of course, be expected to result in the creation of particles when a(t) is changing. But 
here we are only considering the effects of the gravitational metric and its derivatives 
that appear in the generally- covariant free field equation in (I27p . Notice that the 
gravitational field does not need to be weak for the present considerations to apply. 
Therefore, one would expect our results to apply in the early universe for scales below 
the Planck scale. In particular, for the minimally coupled field equation, fl2]), there 
would be strong creation of minimally-coupled scalar particles by the gravitational field 
in the early universe. As noted in [151 US] > this creation process in the early universe 
would react back on the gravitational field through the non-linear gravitational field 
equations. Many studies of this process and its reaction back on the gravitational field 
have since been carried out. 

7. The connection with conformal invariance 

About a year after finding that there is no creation of particles of the field satisfying 
( 1271) for any smooth expansion of the universe, I happened to read an article by Roger 
Penrose [10] , in which he considered conformal transformations of the metric of the form 
g^ v (x, t) — > Q(x, t) 2 g^(x, t) = g^ix-- an d wrote down a field equation of a scalar field 
4>{x,t) that (with multiplication of the scalar field by a suitable power of Q(x, t) to 
define the conformally transformed scalar field 0) remains the same in terms of the 
conformally transformed metric g and field 0. What Penrose stated was that under 
such a conformal transformation of the metric and field the scalar wave equation 

(□+(1/6)^)0 = 0, (31) 

implies that satisfies 

(□ + (l/6)i?)0 = O, (32) 

where □ and R correspond to the □ operator and scalar curvature R, but now calculated 
using the metric g^ v in place of g^ v . 

When I saw that the scalar field equation, (127]) . that gives exactly no creation 
of real particles for any smoothly changing spatially flat FLRW universe was the 
conformally invariant scalar wave equation, I set out to find out why conformal symmetry 
is responsible for the existence of the preferred vacuum state in this class of FLRW 
universes. The reason is the global conformal flatness of the the spatially-flat FLRW 
universes. 

The conformal transformation of the metric that takes one from the spatially flat 
FLRW to Minkowski space is simply given by 

-> a(t)~ 2 v(x,t) = v(x,t), (33) 

followed by defining the new Minkowski-time variable tu = /' a{t')~ 2 dt' . (This 
Minkowski time variable tu is the same as the conformal time variable 77.) It immediately 
follows that for non-singular smooth ait) the conformally transformed equation f l32l 
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with time tM is simply the free field equation in Minkowski space. Hence, because 
there is no creation of real particles in Minkowski space there is none in spatially- 
flat FLRW universes having non-singular smooth scale factors a(t). The positive and 
negative frequency solutions in Minkowski space, when tu is written in terms of t and 
the solutions are multiplied by the appropriate power of the conformal factor that takes 
the Minkowski space <ft t° the FLRW <ft, one finds that the exact positive and negative 
frequency solutions are identical to the positive and negative frequency solutions in fl30|) . 
The existence of a global preferred vacuum state in these FLRW universes, thus follows 
simply from the existence of a preferred vacuum state in Minkowski space. For field 
equations having non-zero m 2 or that break conformal invariance in other ways, there 
is indeed particle creation by spatially-flat FLRW universes. For further details and 
references, see [14J [Chapter 2]. 

8. Adiabatic Regularization of the particle number operator 

Let us return now to the minimally- coupled field of equation (j2J). When I attempted 
to use the particle number operators at(t)a^(t), with a^(t) defined in ffTTj) . to calculate 
the total number of created particles present in a given co-moving volume at time t, the 
following problem arose. Suppose that the scale factor a(t) is a smooth function, with 
[d(t)/a(t)} 2 and/or d(t)/a(t) non-zero. I found that the total density of particles created 
from the initial Minkowski vacuum state was divergent because the number of particles 
created in the co-moving volume did not go to sufficiently rapidly as k approached 
infinity. 

In the continuum limit of large coordinate length L, the expectation value at time 
t, of the number density of particles that were created from the early-time Minkowski 
vacuum state is 



The spatially-flat FLRW universes that were most commonly considered at the time 
were the matter-dominated, radiation-dominated, and steady-state de Sitter universes, 
for which a(t) was proportional, respectively, to t 2 ^ 3 , t 1//2 , and exp(Ht). For these and 
similar expanding universes, dimensionless quantities such as [d(t) / a(t)] 2n [k/a(t)]~ 2n 
and [d(t)/a(t)] n [k/a(t)]~ 2n , for positive integers n, are roughly of the same order-of- 
magnitude, and become small for momenta [k/a(t)} that are large with respect to the 
Hubble parameter, a(t)/a(t). 

Therefore, I expanded the dimensionless quantity 0k(t) in a series of terms having 
increasing numbers of time derivatives of a(t). The number of time derivatives of 
a(t) that appear in a term is called the adiabatic order of the term because terms of 
successively higher adiabatic order approach as higher powers of the Hubble parameter, 
and hence faster in the adiabatic limit of an infinitely slow expansion of the universe. 




(34) 



Here, |0) is the state annihilated by the operators for all modes k. 
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The leading term in the adiabatic expansion of the dimensionless quantity, \(3 k (t)\ 2 , 
turned out to be of second adiabatic order, having dimensionless terms proportional 
to [a{t)/a{t)] 2 [k/a{t)}- 2 and [a{t)/a{t)\ \k/a{t)}- 2 . Therefore, the above integral, (131), 
diverges linearly as the upper limit of integration approaches infinity. This UV 
divergence is not present when the scale factor becomes constant in the late-time 
Minkowski space. It is caused by virtual particles (discussed in the previous Section) 
that are present only while the scale factor a(t) of the universe is changing. 

To deal with this infinity, I considered the physical requirements that should be 
satisfied by a number operator corresponding to the number of real particles that would 
be measured during a time t when a(t) is smoothly changing. As for any ideal mea- 
surement in quantum theory, there must be a corresponding hermitian operator that has 
eigenvalues that correspond to the possible results of the measurement. Because the re- 
sult of an ideal measurement of particle number present at a given time should be an 
integer, I required: 

(a) that the eigenvalues of the hermitian particle number operator, for each indi- 
vidual mode of the field, is the set of non-negative integers; and 

(b) that this hermitian particle number operator would reduce to the number op- 
erator at(t)a^(i) during any periods of time in which the smooth function a(t) became 
constant ( as for example, at early and late times in an asymptotically static expanding 
universe). 

There is an irreducible uncertainty in measuring the particle number present at a 
given time t in a given co-moving volume [a(t)L] 3 in the expanding universe during a 
period when a(t) is smoothly changing and the derivatives of a(t) are non-zero. For 
example, suppose one were to prevent particles from escaping or entering through the 
boundaries of the co-moving volume for a time interval, At, during which one counted 
the number present. If there were no particle creation by the expansion of the universe 
in the co-moving volume, then (as a result of the time-energy uncertainty relation) 
the number counted would have an uncertainty inversely proportional to the time of 
measurement, At. On the other hand, if the particles were being created in the co- 
moving volume by the expansion of the universe during the time At, there would be an 
additional uncertainty directly proportional to At in the number counted. Hence, there 
would be a minimum or irreducible uncertainty in the number of particles counted by 
the observer. 

As we shall see, the irreducible uncertainty in the measured particle number of a 
given quantized field, such as the free scalar field, allows us to define a suitable hermitian 
particle number operator that satisfies requirements (a) and (b) and that has no UV 
divergence. 

Within the theory being considered here, although the linearized perturbations of 
the gravitational field could be treated as a quantum field satisfying a linear differential 
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equation analogous to that of the quantized scalar field under consideration, the 
background FLRW universe itself is regarded as classical. Because of the incompleteness 
of the semi-classical theory, it is perhaps no longer surprising (as it was to me in 1964, 
when I did this work in my Ph.D. thesis) that a regularization and renormalization 
procedure, such as the one to be described below, is necessary to obtain a suitable 
number operator satisfying requirements (a) and (b) for this free (non-interacting) scalar 
field. 

I implemented the requirements (a) and (b) by starting from the number operator 
at(t)a^(t) and defining new creation and annihilation operators by adding terms to 

at(t) and to a^(t) that involved [a(t)] 2 and a(t) in such a way that the new operators 
satisfied, up to second adiabatic order, the commutation relations for creation and 
annihilation operators. This implied that up to second adiabatic order, the new 
hermitian particle number operator formed from these new creation and annihilation 
operators had eigenvalues consisting of the non-negative integers, thus satisfying the first 
requirement. Furthermore, the canonical commutation relations for the quantized scalar 
field and its conjugate momentum were maintained. This was because the definitions of 
the new creation and annihilation operators resulted simply in a renormalization of the 
original frequency u> of the mode k of the field to a new frequency, shifted relative to u 
by a term involving [a(t)] 2 and d(t). This frequency shift is of second adiabatic order. 

These terms of second adiabatic order can be said to arise from virtual particles 
because they do not depend on the past behavior of a(t) and they vanish at early and 
late times in a statically bounded expansion of the universe. These virtual particles 
have the effect of renormalizing the frequencies of the modes of the field. 

In my Ph.D. thesis I obtained expansions to all adiabatic orders of the coefficients 
a.k{t) and (3k(t) in ( fTTl) . These coefficients define a^(t) as a linear combination of the 
constant annihilation and creation operators and Al, respectively. Using these 
adiabatic expansions, I found that I could subtract the terms of second adiabatic order 
from the exact coefficients «&(£) and ^(t) in such a way that the new adiabatically 
regularized operators replacing a^{t) and a>?(t) satisfied the required commutation 
relations for annihilation and creation operators. Hence the eigenvalues of the 
hermitian number operator formed from these adiabatically regularized operators 
satisfied requirement (a) that the eigenvalues of the hermitian number operator be 
the non-negative integers. Requirement (b) is also satisfied because the subtractions 
involve time-derivatives of a{t) and hence vanish for any period during which a(t) is 
constant. Hence, the adiabatically regularized operators do reduce to the original a^(t) 
and al(t) if a(t) becomes constant for a period of time. I also showed that this could be 
continued to higher adiabatic order, but because the expression for the particle number 
had no UV divergence when the second adiabatic order subtractions were made, I felt 
that subtractions of higher adiabatic order were not necessary when regularizing the 
particle number. 

I further showed that the result for the adiabatic subtractions was unique to second 
adiabatic order. For example, one could replace Uk(t), throughout our considerations, 
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by the quantity ^k/a(t)] 2 + m 2 + F(k, a, d, t), where F is a function that vanishes when 
a and a vanish and which never causes that quantity to vanish or become imaginary. 
One would obtain the same second adiabatic-order subtractions required for the creation 
and annihilation operators to satisfy the appropriate commutation relations to second 
adiabatic order. Similarly, to second adiabatic order, one would obtain the same 
result for the regularized quantized field and the renormalized frequency [6J [pp. 63-69, 
Eqs.(80)-(89)]. 

Here I will write down and then confirm the result I obtained for the renormalized 
field from the requirements (a) and (b) by using the adiabatic series expansions for 
a and and showing that the second adiabatic order changes in the creation and 
annihilation operators had the effect of renormalizing the frequency of each mode of the 
field. I will write the end-result for the minimally-coupled scalar field, in terms of the 
regularized creation and annihilation operators and the renormalized frequency for each 
mode. As mentioned above, this procedure gives a particle number operator that has no 
UV infinities when summed over all modes. I will write down the final result below and 
then use a relatively simple method (also given in [6][pp 155-156]) to confirm the result 
to second adiabatic order. This method can be extended to higher adiabatic orders if 
required to renormalize a physical quantity having higher order UV divergences. 

To second adiabatic order, I found that the frequency of each mode is shifted from 
oo to u — S, where S was defined in ( 12X14241) . To confirm that the frequency of each mode 
k is shifted, to second adiabatic order, by S, let us now return to the minimally- coupled 
field equation (j2J) for a particle of mass m. In a spatially-flat FLRW universe, with 
metric 01) , the field equation takes the form (J3]). The functions h k (t) are defined in 
terms of the mode functions fa(x,t) in j5]). The h k {t) satisfy (jHJ), which is written here 
again for the convenience of the reader: 

^h k + ^h k {t) + m 2 h k {t) - | h k {t) - ~h k (t) = 0. (35) 

When we were arriving at the equation for which there is no particle creation (the 
massless conformally invariant equation), we rewrote the above minimally-coupled scalar 
field equation in the equivalent form of (1211) . which we also rewrite here: 

h k (t) + {u k {tf - u k {tf' 2 [{d 2 /dt 2 )uj k {t)- 1 / 2 }} h k (t) = 2u k (t)S(k,t)h k (t), (36) 

where oo k (t) is given in ([TO]) and the function S(k,t) is defined by 

2u k (t)S(k,t) = C x {k,t){a{t)/a{t)f + C 2 (k,t)(a(t)/a(t)), (37) 

with 

fc 4 + 3m 2 a(t) 2 k 2 + (3/4)m 4 a(t) 4 

Cl(M) " (P + mMt) 2 ) 2 (38) 

and 

k 2 + (3/2Ka(t) 2 r3g . 
The function S is clearly of second adiabatic order. 



Particle creation and particle number in an expanding universe 



18 



It is not difficult to confirm that (1361) can be written to second adiabatic order, as 

k(t) + {W k {tf - W^tY^/dt^Wkit)- 1 / 2 }} h k (t) = 0, (40) 

where Wk(t) is here defined as 

W k (t)=u k (t)-S(k,t). (41) 

with u>k(t) defined in ( fTUj) as \J[k/a{t)} 2 + m 2 . The reader will recall from ( 120]) that the 
exact general solution of this equation is ( !T9|) . where ot k and f5 k are complex constants, 
with | ctfc | 2 — \(3k\ 2 = 1. The frequency of mode k to second adiabatic order is thus 
given by W k (t) of (j4~T|) . For this solution having the exact adiabatic form, there is no 
creation of real particles to second adiabatic order because the coefficients ot k and f3 k 
are constants to that order. One can continue to iterate this process to arbitrarily high 
adiabatic order if a(t) has time derivatives to all orders. This does not mean that \f3 k \ 2 
must actually be 0, but only that it must have an essential singularity as a(t) smoothly 
approaches a constant value. An example of a \(3k\ 2 with that property is the function 
exp[— (k / a(t)) 2 R(t)~ 2 }, which has an essential singularity as the time-derivatives of a(t) 
approach 0. Here R(t) is the Ricci scalar curvature of (|26|) . 



Heuristically, this shift in frequency is analogous to the shift in the frequency of a 
pendulum that is undergoing simple harmonic oscillations and is suspended in a viscous 
fluid. The viscosity of the fluid causes a shift in frequency that is analogous to the 
shift in frequency of each mode of the field caused by the presence of short-lived virtual 
particles. Continuing the analogy, if the pendulum bob is suspended from a string, then 
gradually changing the length of the string will tend to excite the pendulum from its 
ground state to higher quantum states, with various probabilities. Such excitations are 
analogous to the creation of particles caused by the expansion of the universe. 

In summary, the minimally-coupled quantized scalar field with frequency 
renormalized to second adiabatic order and the corresponding creation and annihilation 
operators a>t(t) and a k (t) that give a particle number operator N^(t) = at(t)a^(t) that 
has no UV divergence when summed over all modes, is 



where g k (x, t) is 

^(x,t) = K- 1 /V^(2a(t) 3 iy fc (t))- 1 / 2 exp(-z / W k {t')dt'). 

and to second adiabatic order 

W k {t)=u k {t)-S{k,t). 
Also to second adiabatic order, 



af M = a k (t)A s + fo(t) m AL 



with 



(43) 



(44) 



(45) 



0, and 



<**(*). 4(*) ( 46 ) 
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Thus, the o,t(t) and a%(t) satisfy the commutation relations of creation and annihilation 
operators to second adiabatic order at each time t. 

Because S(t) is of second adiabatic order, it vanishes in the adiabatic limit of 
an infinitely slow continuous and smooth expansion of a(t) (in which each successively 
higher derivative of a(t) approaches as one higher power of an adiabatic parameter that 
approaches 0). Thus, in the adiabatic limit, the above expressions P3l) . p5l) and 

f|46|) approach (TT5T) — (T1~8T) . implying that the number operator N^(t) = at(t)a^(t), formed 
using the operator a^(t) in ( 1451) . is an adiabatic invariant. What we have gained by 
working to second adiabatic order is an expression for a particle number operator N^(t) 
that makes (1341) finite at all times t when a(t) is smoothly changing, and which reduces 
to the known exact Minkowski space particle number during any time-interval when a(t) 
has a constant value (as for example, at early and late times for a statically-bounded 
or asymptotically-static expansion of the universe). As mentioned earlier, one can 
continue this process of adiabatic regularization to higher adiabatic orders, depending 
on how many continuous derivatives of a(t) exist. Generally, we prefer to make the 
minimum number of subtractions necessary to renormalize a physical quantity such as 
the expectation value of the particle number, in order to avoid introducing unnecessary 
artifacts of regularization into the predicted results of physical measurements. 

The irreducible uncertainty inherent in the physical process of measuring the 
particle number of the minimally-coupled scalar field has been discussed earlier. That 
irreducible uncertainty implies that we need not expect there to be a unique precisely 
defined particle number operator during an expansion of the universe that is creating 
particles. Hence, there is no contradiction in defining the particle number operator for 
that field only to second adiabatic order. In fact, the adiabatic series for |/3| 2 in general 
is asymptotic, but not necessarily convergent, as explicitly shown in Chapter 2 of [T3] . 

9. A few words on fermion and boson fields of higher spin 

In the article of Penrose [10] that I mentioned earlier, in addition to the conformally- 
invariant massless scalar field equation, he wrote the field equations using 2-component 
spinor index notation for massless fields of higher spin in curved space-time. (The 
2-component massless spin- 1/2 field equation is the same as the Pauli equation for 
a massless neutrino and the spin-1 equation is equivalent to the Maxwell equations 
in vacuum.) Penrose pointed out that these field equations of massless fermions and 
bosons, written using 2-component spinor indices, were all invariant under conformal 
transformations of the metric and fields. In my Ph.D. thesis [6], I used those equations, 
and essentially the same method of proof that I had used for the free, massless, 
conformally-invariant scalar field equation, to show that the free massless fermions 
and bosons described by these conformally-invariant field equations were not created 
by smoothly expanding spatially-flat FLRW universes. These results held for photons, 
massless neutrinos, and the other conformally invariant free massless fields satisfying 
the given conformally-invariant linear field equations of higher spin. In particular, this 
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result also holds for the spin-2 conformally invariant massless quanta (i.e. ,conformally- 
invariant gravitons) that were described by the spin-2 field equation in [TO]. 

However, these conformally-invariant gravitons are not the same as the gravitons 
that are obtained from the linearized Einstein gravitational field equations. In 1946, 
E.M. Lifshitz [2T] had already proved that in the Lifshitz gauge each of the 2 components 
of this linearized "Einstein" graviton field satisfies an equation that is the same as 
that of the massless minimally- coupled scalar field equation. I did not know about the 
1946 paper of Lifshitz [21] . so in j6j H5J HE] I reported on my result that minimally 
coupled scalar particles would be created by the FLRW expanding universe, and I also 
reported on my result that free massless fields satisfying conformally-invariant field 
equations would not be created by these FLRW universes. In particular, I reported that 
the spin-2 gravitons that satisfied the conformally invariant massless spin-2 free field 
equation would not be created in such a universe. Not being aware of Lifshitz graviton 
field equation obtained from the Einstein gravitational field equations, I did not use 
my proof for the minimally-coupled scalar field to report that linearized "Einstein" 
gravitons would indeed be created by an expanding FLRW universe. The creation of 
these linearized Einstein gravitons was studied some time later [221 123] . 

As mentioned above, through an argument based on conformal invariance one can 
show that the free massless spin- 1/2 Dirac equation would not give rise to creation 
of massless neutrinos in an FLRW universe, and that the quantized free Maxwell field 
equation does not give rise to the creation of photons from vacuum in an FLRW universe. 
Similarly, the presence of free photons (that originated from distant sources such as 
stars) in intergalactic space would not give rise to stimulated creation of photon pairs 
simply as a result of the expansion of an FLRW universe 0, [15j [16], [24] . 

This brings us to some remarks on work done by Schrodinger in 1939 and 1940 
[TT| [12]. He considered a quantum mechanical wave packet that is propagating in an 
expanding FLRW universe. He showed that such a wave packet would increase in 
amplitude and would also produce a weaker backward-going wave packet as a result of 
the expansion of the universe. He interpreted this process as implying that a scalar 
particle (which he considered in his actual derivation) or, more importantly, a photon 
moving in a given direction had a non-zero probability of inducing the creation of 
a pair of photons, one moving in the same direction as the original photon (which 
continues propagating in its original direction) and the other moving in the opposite 
direction. The frequency of the created photons would be the same as that of the original 
propagating photon. This process of pair creation is stimulated by the existence of a 
particle propagating in a given direction. In the case of photons, he suggested that light 
propagating from a source, such as a distant star, would create more photons moving in 
the original direction (as well as the others member of each created pair moving in the 
opposite direction). At the time when the pair is created, the created photons would 
have the same frequency as the original photon, and that frequency would be red-shifted 
by the expansion of the universe, just as the original photon's frequency is shifted in 
the FLRW expanding universe. He called this an "alarming" phenomenon. 
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In the quantum mechanical context in which Schrodinger was working, this process 
is the stimulated creation of pairs of photons from the presence of existing photons in 
the expanding FLRW universe. However, as mentioned earlier in this article, I showed 
that because photons in the FLRW expanding universe satisfy a conformally invariant 
equation, there is no creation of photon pairs by the expansion, either of the stimulated 
variety that he was considering, or of the spontaneously created pairs from the vacuum. 

I had been unaware of this work of Schrodinger on stimulated pair creation until the 
early 1970's, when a former colleague of Schrodinger told me about it. Soon afterward, 
I wrote a paper pointing out what I said above [23]. Thus, Schrodinger need not have 
been alarmed by this process. 

Since there seems to be some confusion about the pair creation process from the 
vacuum, I should point out that the pair creation of scalar particles from vacuum was 
first discovered in my Ph.D. thesis [6]. The method of adiabatic regularization that I 
have discussed in the present paper was developed to deal with UV infinities that arise 
in this process of pair creation from the vacuum during the time when the scale factor 
of the universe, a(t), is actually changing with time. The stimulated creation of pairs is 
also discussed in P, [151 OS] • 

10. Conclusion 

My main purpose in writing this paper is to explain in some detail how the method 
of adiabatic regularization stemmed from the fundamental physical requirements that 
must be met by a particle number operator in an expanding FLRW universe. This 
development was prompted by the discovery [6] that particle pairs were created from 
the vacuum in an FLRW universe with a smoothly changing scale factor a(t). The 
number of particles created in any co-moving volume was finite for each mode of the 
quantized field. Furthermore, the total number summed over all modes was finite if one 
counted them after the universe had smoothly stopped expanding. 

However, at any given time when the universe was actually expanding, there seemed 
to be a UV divergence in the total number density of created particles summed over all 
modes of the field, even if a(t) were slowly changing. As observations show that we do live 
in an expanding universe, it became necessary to reconsider the particle number operator 
during times when a(t) is smoothly changing. As we have seen, adiabatic regularization 
started from the physically motivated assumptions that the particle number of the free 
field in each mode (a) should be an adiabatic invariant, analogous to the quantum 
number of the harmonic oscillator with a smoothly changing frequency, and (b) should 
have the non-negative integers as its spectrum of eigenvalues. The existence of particle 
creation caused by the expansion of the universe implied that the measurement of the 
particle number during the expansion of the universe would have a minimum uncertainty 
that depended on the time-interval over which the measurement was made and on the 
rate at which particles were being created by the expansion of the universe during the 
measurement. This allowed the possibility of defining a satisfactory number operator 
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during the expansion of the universe, up to a given adiabatic order (defined in terms 
of number of derivatives of a(t), or equivalently, in terms of powers of an adiabatic 
parameter that goes to when the expansion parameter approaches a constant). It 
turns out that the requirements above could be met, while keeping the total number 
density of particles finite when summed over all modes. That is the main subject of 
this paper, and is being explained in detail here because it has not been fully discussed 
in the literature, outside of [6]. In the course of the present discussion, I also take 
the opportunity to discuss the related topics of conformal invariance and higher spin 
quantized fields. 

In view of Professor Stuart Dowker's interest in renormalization in curved space- 
time and his outstanding contributions, it seems fitting to present these results in this 
special issue in his honor. 

References 

[1] Dowker, S J and Critchley, R 1976 Phys. Rev. D13, 3224. 

[2] Dowker, S J and Critchley, R 1977 Phys. Rev. D16, 3390. 

[3] Dowker, S J and Banach, R 1978 J. Phys. A: Gen. Phys., 11, 2255. 

[4] Hawking, S W 1977 Commun. Math. Phys. 55, 133. 

[5] Parker, L 1979 "Aspects of Quantum Field Theory in Curved Spacetime: Effective Action and 
Energy-Momentum Tensor," in Recent Developments in Gravitation: Cargese 1978, edited by 
M. Levy and S. Deser (Plenum Publishing Corp., New York), pp. 219-273. 
[6] Parker, L 1966 The Creation of Particles in an Expanding Universe, Harvard University, 

Cambridge, Massachusetts. 
[7] Glenz, M M and Parker, L 2009 Phys. Rev. D 80, 063534. 
[8] Parker, L and Fulling, S A 1974 Phys. Rev. D9, 341. 
[9] Fulling, S A and Parker, L 1974 Ann. Phys. (N.Y.) 87, 176. 
[10] Penrose, R 1964 in Relativity, Groups and Topology, edited by C. DeWitt and B. S. DeWitt 

(Gordon and Breach, New York), 565, 566. 
[11] Schrodinger, E 1939 Physica 6, 899. 
[12] Schrodinger, E 1940 Physica 46, 25. 

[13] Birrell, N D and Davies, P C W 1982 Quantum Fields in Curved Space (Cambridge University 
Press, London) 

[14] Parker, L and Toms, D J 2009 Quantum Field Theory in Curved Spacetime: Quantized Fields and 

Gravity (Cambridge University Press, London) 
[15] Parker, L 1968 Phys. Rev. Lett. 21, 562. 
[16] Parker, L 1969 Phys. Rev. 183, 1057. 
[17] Parker, L 1971 Phys. Rev. D3, 346. 
[18] Agullo, I and Parker, L 2011 Phys. Rev. D83, 063526. 

[19] Agullo, I and Parker, L 2011 Gen. Relativ. Gravit., DOI: 10.1007/sl0714-011-1220-8. First Award 

in the 2011 Essay Competition of the Gravity Research Foundation. 
[20] Agullo, I, Navarro- S alas, J and Parker, L 2012 "Enhanced local-type inflationary trispectrum from 

a non- vacuum initial state" JCAP (accepted for publication). 
[21] Lifshitz, E M 1946 Zh. Eksp. Teor. Fiz. 16, 587. 

[22] Grishchuk, L P 1974 Zh. Eksp. Teor. Fiz. 67, 825; Sov. Phys. JETP 40, 409 (1975). 
[23] Ford, L H and Parker, L 1977 Phys. Rev. D16, 1601. 
[24] Parker, L 1972 Phys. Rev. D5, 2905. 



